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Abstract 

We report on density waves in granular material, investigated both ex- 
perimentally and numerically. When granular material falls through 
a long narrow pipe one observes recurrent clogging. The kinetic en- 
ergy of the falling particles increases up to a characteristic threshold 
corresponding to the onset of recurrent clogging and density waves of 
no definite wavelength. The distances between regions of high density 
depend strongly on the initial conditions. They vary irregularly with- 
out any characteristic time and length scale. The particle-flow was 
investigated using 2D Molecular Dynamics simulations. Experimental 
investigations lead to equivalent results. 



PACS. numbers: 05.60, 47.25, 46.10, 02.60 



The astonishing effects observed in moving granular material have been sub- 
jects of interest for at least 200 years. Examples of such effects are the angle 
of repose p size segregation f^-p, heap formation on vibrating media 



density waves and granular flow down inchned surfaces [|T^]. 
The origin of these effects is that granular material can behave like a solid or 



like a liquid, depending on its density |T3|. Particularly, in recent times dy- 
namic as well as static behaviour of granular material has been investigated 
by many authors experimentally and theoretically with various techniques 
either analytical, such as thermodynamic fl^ and hydrodynamic ap- 



proaches, or numerical like cellular automata [|T6|, Monte-Carlo simulation 
[^] and molecular dynamics |^ p| Ijl^. Even random walk approaches have 
been proposed [jl8| . 

The aim of this Letter is to report the effect that granular material like 
sand that falls through a vertical narrow pipe or capillary with a diameter of 
only few particle-diameters changes from a homogenous to an inhomogenous 
flow when the kinetic energy Ek of the falling grains reaches a characteris- 
tic threshold. In this state one observes recurrent clogging of the particle 
flow (stick slip motion) and density waves. We simulated this behaviour 
using molecular dynamics and demonstrate that the results agree with our 
experimental observations. 

The experiments were carried out in two vertical fixed glass-pipes of di- 
ameters Di = 2 mm and D2 = 4: mm and length L = 1.4 m. An upper 
funnel contained enough granular material to ensure time-independend ini- 
tial conditions. During the experiments very fine-grained sand of typical 
particle diameter Dg = 0.18 mm begins to flow homogeneously with the 
velocity f = at the upper end of the pipe. After a typical distance of ap- 
proximately 20 cm the uniform flow becomes unstable and turns into a flow 
of recurrent clogging. In this regime density waves can be observed. Fig. |I| 
shows the pipes with typical density distributions. The distances between 
dark high-density regions vary irregularly. There is no definite wavelength. 
The distance between the top of the tubes and the onset of the stick slip mo- 
tion as well as well as for the distances between the plugs depends sensitively 
on the humidity of the air and varies significantly from one experiment to the 
next one, i.e. depend on initial conditions. According to our observations 
the distribution of the distances appears to obey a power law. 

In fact we have done the simulations first and the reported experiments 
were intended to check the accuracy of the results predicted by the numerical 



simulations. 

According to the experiments reported above we investigated numerically 
the flow of spherical-shaped grains with rotational degrees of freedom falling 
under gravity through a pipe of diameter D (fig. ^ using 2D molecular 
dynamics simulations. The length of the pipe considered was much larger 
then its diameter L ^ D. We investigated systems of particles of equal radii 
Ti = R as well as of particles with a Gaussian probability distribution for the 
radii with mean value R (fig. For the flow we assumed periodic boundary 
conditions. The inner surface of the pipe is build of smaller spheres with a 
diameter = | ■ i? of the same material as the grains to simulate a rough 
surface. Since our model includes elastic and inelastic interaction between 
the grains each other and between the grains and the wall we simulated the 
wall exactly in the same manner as the grains but the particles the pipe is 
build of were neither allowed to move nor to rotate. 

The grains with mass rrii were accelerated due to gravity {g = —9.81 m/s"^, 
their density ^^'3 is set to unity): 

= -rrti ■ g (1) 

where is the force acting upon the ith particle in the direction of the axis 
of the pipe. 

Starting with given initial positions of the particles and randomly chosen 
small initial linear and angular velocities the dynamics of the system was 
determined by integrating Newtons equation of motion numerically for each 
particle for all further time steps. For the integration we used a sixth order 
predictor-corrector method for the particle positions and a fourth order one 
for the calculation of the particles' angular movement . 



In our model particles are allowed to penetrate each other slightly. The 
force acting between two particles i and j is assumed to be zero in case 
the particles do not touch each other, i.e. when the distance between their 
centres is larger then the sum of the particle radii Vi + rj. Otherwise the force 
separates into normal and shear force (eq. 0). 



otherwise 



with 



Fn = kN ■ {Vi + Tj - \Xi - Xj\y-^ + 7Ar ■ TUeffiXi - Xj) (3) 

and 



where 



Fs = min{-75 ■ nieff ■ Vrei , ■ \Fn\} (4) 



Vrel = {Xi- Xj) + R - {Qi - Qj) (5) 



rrii ■ rrij m , , 

m^(( = — = — (6) 



The terms x,, Xi and Qi stand for the coordinate, the velocity and the angular 
velocity of the zth particle. Eq. |^ includes the Hertzian contact force which 
rises with the power 1.5 of the penetration depth rj + I of two 



particles. It is used to mimic 3D behaviour of the grains pO|. This ansatz 



for the force was suggested in |^ and slightly modificated in 0. Eq. ^ 



takes into account, that the maximum momentum two particles are able to 
transmit while colliding is determined by the Coulomb-law |]T[. /i is the 
Coulomb-coefficient of the particles. 

To discuss the observed phenomena we refer to simulations with = 600 
particles. The length of the pipe was L = 666 -i? and its width D = 5-R. The 
shear friction coefficient was 7s = 3 ■ 10^s~^, the normal friction coefficient 
7Ar = 3 ■ 10^s~^, the material constant = IO^j^tf and /i = 0.5. For the 
integration time step was chosen At = 10~^ s. 

After filling the particles into the pipe they start moving at rest. While 
falling they gain kinetic energy due to the constant acceleration. To avoid 
infinitely rising velocity of the particles without colliding they were initialized 
with a very small random velocity perpendicular to the axis of the pipe. At a 
certain energy the fiow does not get faster anymore but the grains are decel- 
erated very intensively. From now on the system moves with approximately 
constant kinetic energy. Fig. ^ shows the evolution of the kinetic energy of 
the transversal particle movement Ej. and the energy of the particle rotation 
Ej.. 

In the case of the low velocity regime the granular material moves through 
the pipe with almost homogeneous particle density. Fig. ^ displays snapshots 
of the pipe each 1,000 time steps. Time increases upwards while horizontally 



one sees the evolution of the density wave from the left to the right. Gravity 
acts from left to right. In the early times the particle-density is approximately 
homogeneous due to the initial conditions. During the evolution it gets more 
and more inhomogeneous and density waves are observed. As visible at later 
times the flow is unstable: the regions of high density can diverge as well 
as converge while the average velocity remains approximately the same as 
shown in fig. ^. 

After the transition into the recurrent-plug-flow regime, where the ve- 
locity profile does not vary with the distance from the centre on the pipe, 
the system reaches a dynamic state of nearly constant kinetic energy due 
to the equilibrium between constant acceleration and energy dissipation by 
friction. Nevertheless one can observe small fluctuations of the energy which 
are aproximately due to the genesis and vanishing of regions of higher density 
with time, as visible comparing fig. ^ with fig. ^. This behaviour proves that 
the inhomogenous flow is unstable and it could be an indication of coexisting 
metastable states but we did not check this possibility yet. Our simulations 
show further that there are no significant differences between the behaviour 
of a system with particles of equal and randomly distributed radii. Simu- 
lating the system with slightly different initial conditions we got in principle 
the same results for the energy evolution and for the structure of the density 
wave but the concrete density wave as it can be observed as dark and bright 
regions in fig. ^ varies significantly. Very small differences in the initial con- 
ditions lead to very different macroscopic behaviour. In this sense we call 
the flow chaotically. 

There are essential differences between the behaviour of a liquid and a 
fluidized granular material: for the case of an incompressible fluid one never 
flnds density waves. Nevertheless there are similarities too: The evolution 
of the energy of the system is similar to the evolution of the velocity of 
a liquid flowing through a pipe [^] while the pressure is increased gently. 
At low flow velocity an incompressible viscous fluid forms a laminar flow 
(Poisseuille-flow) which gets faster with rising energy. If the velocity reaches 
a certain value Uc the Poiseuille-flow becomes unstable and an inhomoge- 
neous flow regime becomes stable. For the case of a flxed pipe geometry and 
given fluid parameters the critial velocity of the flow Uc is described by the 
critical Reynolds-number p3| Rcc = ~ 2, 300 where Hk is the kinematic 
viscosity of the fluid. Provided there are critical fluctuations the system turns 
into the turbulent regime. If the pressure is increased very gently and the 



system is not disturbed in another way, however, it is possible to maintain 
the Poiseuille-flow at velocities far over Rcc in the unstable regime. Then a 
small fluctuation causes the system to turn to the stable state, i.e. it sud- 
denly lowers the velocity of the flow and dissipates the energy. As shown for 
the flow of a granular material there exist different flow regimes too and the 
transition between them is due to a sudden energy dissipation. Hence the 
transition from low energy to high energy regime of granular flow in a pipe 
is similar to the transition of a laminar fluid-flow into a turbulent one. 

Experimental and recent numerical investigations of traffic flows lead 
to similar results as those presented here. At a given average car-velocity 
the homogeneous traffic flow becomes unstable and traffic jamming occurs 
provided there are random perturbations of the velocity of the cars which 
correspond to the collisions of the grains with the wall in our case. 

Concluding we state that granular material falling through a narrow pipe 
appears effects far from the behaviour of a liquid. At a certain energy the 
granular flow transits from the homogeneous to an inhomogenous flow-regime 
where clogging and density waves can be observed. The regions of high par- 
ticle density can converge as well as diverge, their distances vary irregularly. 
The results of the experiment shown in the beginning agree with the numer- 
ical simulations. 

The author thanks J. Gallas, H. Herrmann, G. Ristow and S. Sokolowski for 
fruitful and stimulating discussions and H. Puhl for experimental assistance. 
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List of Figures 

1 The experiment: fine grained sand with a typical particle di- 

ameter of Dg = 0.18 mm flows through a glass capillary tube. 

Recurrent clogging and density waves as well as free fallinc, 
grains are visible. The distances between the (dark) regiom 
of high density fluctuates significantly. The figure shows two 
pipes with different diameters [2 mm and 4 mm, length I.4 m\ 

and an enlargement of the thinner one. 

p The pipe. The rough surface consists of spheres with diameter 

Tg = \ ■ H of the same material as the grains. 

3 Enlargement of a small part of the pipe filled with grains of 

different radii. Gravity acts horizontally from left to right. . 

1 Evolution of the kinetic energy of the transversal particle move- 
ment and the rotational energy. Since the rotational energy it 
very small the curve is almost indistinguishable of the time- 
axis. While ffowmg accelerated, the system becomes unstable 
at 35.000 time steps approximately and transits into another 

flow regime where it reaches a state of relative steady kinetic 

energy (E ^ 5 ■ 10'^ for this simulation) .\ 

f) Snapshots of the pipe each 1,000 time-steps. Time increase 

upwards. The density fluctuations are not equidistant and not 

stable. 'I'hey can converge as well as diverge. Dark regiom 



correspond to higher densities 



Figure 1: The experiment: fine grained sand with a typical particle diameter 
of Dg — 0.18 mm flows through a glass capillary tube. Recurrent clogging and 
density waves as well as free falling grains are visible. The distances between 
the (dark) regions of high density fluctuates significantly. The figure shows 
two pipes with different diameters (2 mm and 4 mm, length I.4 m) and an 
enlargement of the thinner one. 
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Figure 2: The pipe. The rough surface consists of spheres with diameter 
rs — \ ■ R of the same material as the grains. 



Figure 3: Enlargement of a small part of the pipe filled with grains of different 
radii. Gravity acts horizontally from left to right. 



Figure 4: Evolution of the kinetic energy of the transversal particle movement 
and the rotational energy. Since the rotational energy is very small the curve 
is almost indistinguishable of the time-axis. While flowing accelerated, the 
system becomes unstable at 35. 000 time steps approximately and transits into 
another flow regime where it reaches a state of relative steady kinetic energy 
(E 5 ■ 10^ for this simulation). 



Figure 5: Snapshots of the pipe each 1,000 time-steps. Time increases up- 
wards. The density fluctuations are not equidistant and not stable. They can 
converge as well as diverge. Dark regions correspond to higher densities. 



